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We have simulated, using parallel tempering, the three dimensional Ising spin glass 
model with binary couplings in a helicoidal geometry. The largest lattice (L = 20) has 
been studied using a dedicated computer (the SUE machine). We have obtained, measuring 
the correlation length in the critical region, a strong evidence for a second-order finite 
temperature phase transition ruling out other possible scenarios like a Kosterlitz-Thouless 
phase transition. Precise values for the v and n critical exponents are also presented. 



I. INTRODUCTION 

The study of spin-glasses [lj, beyond their own physical relevance, has opened new ways to Statistical 
Physics. The solution Q of the Sherrington-Kirkpatrick model, that describes spin-glasses living in infinite 
dimensions, allowed the introduction of a new set of ideas that have found applications in very different con- 
texts, like Optimization, Neural Networks, and so on. Yet the applicability of the rich infinite-dimensional 
physical picture to describe the low-temperature physics of three-dimensional spin-glass materials (like, for 
instance, CuMn, AgMn, Eu^ Sri^ S, see S) is still controversial [Q. Furthermore, a rather simpler ques- 
tion, what is the nature of the spin-glass phase transition?, has not yet found a fully satisfactory answer. 
Although the very existence of a phase transition has been questioned, from the experimental side, there 
is now a wide consensus on its existence, as signalled by the behavior of the non- linear susceptibility |5). 

On the other hand, the theoretical approach is almost limited to the Monte Carlo simulation of the 
Edwards- Anderson model, given the enormous difficulties found when using field-theoretic renormalization 
group techniques (|]. Recent numerical simulations RjO] have found indications of a finite-temperature 
phase transition, which has been confirmed in Ref. Jl2|T However, the possibility of a Kosterlitz- 
Thouless like phase transition (an exponential divergence of the correlation-length at the critical temper- 
ature followed by a line of critical points) could not be excluded p 2 | . Even so, critical-exponent estimates 
that could be compared with experiments were obtained [[7H15[by assuming power-law divergences at 
the critical temperature (i.e., non Kosterlitz-Thouless behavior). However, the statistical errors of those 
estimates (10% for the correlation-length exponent v, and 15 % in the anomalous dimension, rj), and that 
of the critical temperature estimate seems poor compared to similar computations for ordered systems, 
which is due to the numerical difficulties encountered on the simulation of the Edwards- Anderson model. 
In fact, the issue (crucial for accurate calculations of critical exponents) of the scaling corrections has not 
been addressed in previous works, exception made of Ref. [ jl2| . 

In this work, we shall perform a detailed study of the critical behavior of the Edwards- Anderson model. 
The numerical simulations have been in part performed on a dedicated computer (the SUE machine, see 
below for more details), on which we have been able to thermalize 6920 samples of 20 3 lattices at the 
critical temperature, the largest thermalized lattices in previous studies at similar temperatures being 16 3 . 
For the thermalization deep inside the critical region we have used the Monte Carlo exchange method (also 
known as parallel tempering) |lq-[lq| . 

Our study shares with Ref. |12||~the definition of the finite-lattice correlation-length |]l9| , and a heavy 
use of the Finite-Size Scaling (FSS) Ansatz p3). Yet, both analyses are rather different. Ref. [O] uses 



the techniques of Ref. |21| to extrapolate the measures taken on lattices which are small compared with 
the correlation- length, to the thermodynamic limit. On the other hand, we use the quotient method p2| , 
where measures taken on two lattices are compared at the temperature at which the correlation-length 
measured in units of the lattice-size coincides for both. 



For the particular problem of spin-glasses, the method of Ref. has the advantage of not requiring the 
thermalization of large lattices at large correlation-lengths (the dynamical critical exponent for the three 
dimensional Ising spin glass in the critical region is near 7 p3[ ). On the other hand, the quotient method 
offers the possibility of extremely precise determinations of critical exponents and temperatures, and a 
rather transparent control of scaling-corrections, also in disordered systems p4| . The main drawback for 
its use on spin-glass systems is that it requires measures taken on several pairs of lattice of widely different 
sizes at the critical temperature, which is rather difficult due to the above mentioned thermalization 
problems (25). 

We obtain very precise estimates for critical exponents which are compared with the estimates of other 
groups and with available experimental results. The issue of scaling-corrections will be discussed, and a 
rough estimate of to (the correction-to-the-scaling exponent) will be obtained. 

An additional bonus of our computational strategy is that high-quality data for the spin-glass correlation 
length are generated on large lattices at the critical region. This allows for a detailed comparison with 
the archetypical model displaying a Kosterlitz-Thouless phase transition, the XY model in two dimensions 
(for which a cluster method is available, making the simulation almost costless). The finite-size scaling 
behavior of both models is not only quantitatively but qualitatively different. We therefore present strong 
evidence against Kosterlitz-Thouless behavior on the Edwards- Anderson model. 

Finally, we also consider the question of the appropriate cumulant for the study of the spin-glass phase 
transition. In Ref. |p6| it has been argued that the Binder cumulant 27 works poorly, in marked contrast 
with ordered systems. It is also claimed that the cumulant G introduced in Ref. |2§| for the study of 
systems without time-reversal symmetry, see Eq.(Q), does a better job. We shall show that it suffers from 
similar scaling corrections but of opposite sign, so that its crossings happens at temperatures higher than 
the critical point. This is rather advantageous from the point of view of thermalization. On the other 
hand, its measures are far noisier than the ones of the finite-lattice correlation-length, and it also suffers 
from stronger corrections to scaling. 

The large statistics needed to obtain precise results on larger lattices has been possible by the use of 
a dedicated computer based on programmable components. Details about the machine can be seen in 
Refs. |3f|. SUE, for Spin Update Engine, consists of twelve boards attached to a PC. Each of these 
boards contains programmable devices and memories, allowing the simulation of eight lattices of size up to 
L = 60. The presence in each board of a component dedicated to random number generation allows the use 
of the Heat Bath algorithm. During the simulation, Schwinger-Dyson equations (3^]3^] were used to check 
for the correctness of the random-number sequences. Periodically, the spin configurations are downloaded 
to the PC, where the measures are performed and stored. When using the parallel tempering scheme, 
the PC controls the mechanism, interchanging the configurations corresponding to adjacent temperatures 
when appropriate. The update speed of the whole system is 0.22 ns/spin, one hundred times faster than 
one Alpha EV5, 600 MHz processor running multi-spin code. 

The layout of the rest of this paper is as follows. In the next section we introduce the model, the definition 
of the observables and we review the basis of our Finite Size Scaling (FSS) method. The third section 
is devoted to analyze the statistical quality of our data, in particular we discuss the thermalization, the 
parallel tempering parameters and our choice of the number of samples versus number of steps inside one 
sample. The discussion of our numerical results follows, where we compare them with previous numerical 
simulations and experiments. We end the paper with the conclusions. 



II. MODEL, OBSERVABLES AND THE FSS METHOD 

We have studied the three dimensional Ising spin glass defined on a cubic lattice (L x L x L) with 
helicoidal boundary conditions |34j , whose Hamiltonian is 

H = - ^ CTiJijCTj . (I ) 

<i,j> 

The volume of the system is V = L 3 , Oi are Ising variables, Jij (uncorrelated quenched disorder) are ±1 
with equal probability, and the sum is extended to all pairs of nearest neighbors. Due to the quenched 
nature of the disorder, one needs to perform first the thermal average for a given configuration of the Jy 
(denoted by (• ■ ■) hereafter), and later the average over the disorder realization (that will be indicated by 



2 



an ovcrlinc). The choice of hclicoidal boundary conditions is mandatory (for us) because the hardware of 
the SUE machine has been optimized for them. 

We have simulated the smaller lattice sizes (L = 5 and 10) in parallel machines built of Pentium-Pro 
processors (the RTNN machines) using multi-spin coding. We have checked that the L = 10 and L = 5 
lattices are properly thermalized with a standard Heat Bath method (without parallel tempering). The 
larger lattice [L = 20) has been simulated in the SUE machine using parallel tempering and Heat Bath. 

We shall describe in depth the thermalization test and the total statistics achieved in the next section. 



A. Observables 



It is well known that observables in spin-glasses need to be defined in terms of real replicas, that is, 
for every disorder realization, one considers two thermally independent copies of the system {s^, s^}. 
Observables are most easily defined in terms of a spin-like field, the so-called overlap field: 

* = W>. (2) 
The total overlap is the lattice average of the q% 

i 

while the (non-connected) spin-glass susceptibility is 

X q = VW)- (4) 

In Finite-Size Scaling studies, it is useful to have dimensionless quantities, that go to a constant value at 
the critical temperature. A standard example of such a quantity is the Binder cumulant 

Other example is the g2 cumulant [Q , that measures the lack of self-averageness of the spin-glass suscep- 
tibility 



{<, 2 Y 

Of course, any smooth function of these two dimensionless quantities, g-i and <?4, is dimensionless itself. In 
Rcf. |28| it was proposed to study the cumulant G 

because it exhibits a significant reduction of scaling corrections. It has also been argued in Ref. |p6[ that G 
can be extremely helpful for the characterization of the spin-glass phase, but this issue is out of the scope 
of this work. However all the above defined dimensionless quantities, 54, gi and G, require the evaluation 
of a four-point correlation function, which is statistically a much noisier quantity than a two-point one. 
One observable of this kind is the correlation-length, which is defined in terms of the two-point correlation 
function, and its quotient with the lattice size is again dimensionless. Wc therefore arc faced with the 
problem of defining a correlation-length on a finite lattice. In Ref. Jig) ] it was shown how to do it, from 
the following considerations. Let us call C(r) the correlation function of the overlap field, 

CW = ^M (8) 

i 

and C(k) its Fourier transform. Notice that the spin glass susceptibility is simply C(0). Then, inside the 
critical region on the paramagnetic side and in the thermodynamical limit, one has 
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C(k) oc 



1 



2 = _]_dC_ 
C dk 2 



Wl«r 1 



(9) 
(10) 



fe 2 =0 



On a finite lattice, the momentum is discretized, and one uses jl9| a finite-differences approximation to 
Eq.©, 



4 [sin 2 (fc* /2) + sin 2 (fc^/2) + sin 2 (fc^/2)] 



Xq 



C(k m ) 



(11) 



where Xq was defined in Eq. (^) and fc m is the minimum wave- vector allowed for the used boundary con- 
ditions, which in our case is fc m = (2n/L, 2n/L 2 , 2ir/L 3 ). Of course, Eq.jicj) holds on the thermodynamic 
limit (L ^S> £) of the paramagnetic phase. As we do not use connected correlation functions, £ has sense 
as a correlation length only for T > T c . 

We can study the scaling behavior of the finite-lattice definition ( |i~l| ) on a critical point, where the 
correlation function decays (in D dimensions) as r~( D ~ 2+v \ The behavior of the Fourier transform of the 
correlation function for large L in three dimensions is given by 



C(k) 



dr t 



_ 7) sin(fcr) 
kr 



(12) 



and one finds that Xq/C(k m ) goes to a constant value, larger than unity, because ||fc m || = 0(1/ L). 
Furthermore, £/L tends to a universal constant at a critical point (like the Binder cumulant g^j. Moreover, 
on a broken-symmetry phase, where the fluctuations of the order parameter are not critical, one has 
Xq = 0(V), while C(k m ) = 0(1). Therefore the full description of the scaling behavior of £/L is as 
follows. Let £oo be the correlation-length on the infinite lattice: in the paramagnetic phase, for L 3> 
one has £/L = 0(1/ L). On the FSS region, where ^ > L, £/L = 0(1), while on a broken-symmetry 
phase on a lattice larger than the scale of the fluctuations, £/L — 0(L D / 2 ). Consequently, if one plots 
£/L for several lattice sizes as a function of temperature, the different graphics will cross at the critical 
temperature. 

Finally, a very useful quantity is the value of the Hamiltonian, which is used to measure the derivatives 
of a generic observable, O, with respect to the inverse temperature (3 



d p (o) = (on {1) + on (2) ) - (o)(h {1) +n (2) ) 



(13) 



where the W and ( 2 ) superscripts refers to the two replicas needed to construct the operator O. Similarly 
one generalizes the standard reweighting method |55| , that allows to extrapolate the measures taken at (3 
to neighboring values: 



(O) 



/3+A/3 



(O 



^ e A/3W( 1 )+A / 3'H( 2 )^ ' 



(14) 



B. The FSS method 



We have used the quotient method |22|, in order to compute the critical exponents. We recall briefly the 
basis of this method. Let be a quantity diverging in the thermodynamical limit as t~ x ° (t = T/T c — 1 
being the reduced temperature). We can write the dependence of O on L and t in the following way ppf 

0(L,t) =L X °I 

where Fo is a (smooth) scaling function and (— ui) is the corrections-to-scaling exponent (e.g., to is the 
largest negative eigenvalue of the Renormalization Group transformation) . This expression contains the 



Fn 



£(oo,t) 



.0(L- u ,C w ) 



(15) 
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not directly measurable term £(oo, t), but if we have a good definition of the correlation length in a finite 
box £(L,i), Eq. (Il5|) can be transformed into 



0(L,t) = L x °^ 



Co I +OU-) 



(16) 



where Go is a smooth function related with Fo and F% and the term has been neglected because we 
are simulating deep in the scaling region. We consider the quotient of measures taken in lattices L and sL 
at the same temperature 

n / T .s Q(sL,t) 

Qo{s ' L ' t) = ^W (17) 

Then, the main formula of the quotient method is 

Qol Q?=s = s xo/,y + 0(ir w ) , (is) 

i.e., we compute the reduced temperature t, at which the correlation length verifies £(sL, t)/£(L, t) = s 
and then the quotient between 0(sL,t) and 0(L,t). In particular, we apply formula ( |18| ) to the overlap 
susceptibility, \qi an d the- /3-derivative of the correlation length, <9^£, whose associated exponents are: 

Xd p £ = 1- + V, (19) 

x x? = (2 - n)v. (20) 

Notice that Qo\q ( = s can be measured with great precision because of the large statistical correlation 
between Qo and Q^. It is very important that in order to use Eq.([l8|) one does not need the infinite- 
volume extrapolation for the critical temperature but, instead, a reweighting method is crucial to 
fine-tune the = s condition. From Eq.(|lS|) one directly extracts effective exponents (i.e., lattice-size 
dependent) and later on checks for scaling-corrections. Another advantage of the quotient method is that 
the crossing temperature for £/L (i.e., the temperature for which = s) scales as 

grossing K L - u -l/u ^ (21) 

so that one works a factor closer to the critical point than with other FSS methods, such as measuring 
at the maximum of (say) the connected susceptibility. This puts considerably less stress on the quality of 
the reweighting method. It should also be mentioned that one could modify Eq. (|l8[) , and measure at the 
crossing point of a cumulant such as (74, 172 or G. The results should coincide, up to scaling corrections, 
but given the better statistical quality of the measures of £/L, the error bars would significantly grow. 



III. STATISTICAL QUALITY OF THE DATA 

When designing a simulation for a disordered model, one needs to carefully consider how many measures 
will be taken on each sample, N\, and the number of samples to be simulated No,- Two competing effects 
need to be balanced for this. On the first place, from the error analysis of a generic observable, O, one has 
(r) is a Gaussian number of zero-mean and unit variance) 

((0) M ° nteCarl ° - JO)) 2 = ^ (<4o + ^) , (22) 

where as,0 is the variance between different samples of the exact thermal averages, cr^o is the disorder- 
averaged variance for the measures on a sample, and to is an averaged (integrated) autocorrelation 
time [Q. This shows that the optimum value of iVi cannot be much greater than 2too\ j a\. On 
the other hand, when evaluating non-linear functions of thermal-averages, as in Eqs.([l3|,[l4]), a bias of 
order 2to/N\ is present (for the reweighted measures a bias polynomial in 2to/N\ is expected). If \/Ng 
is not much smaller than N\, the statistical errors will shrink enough as to uncover the bias, and we have 
two conflicting goals for the optimization of N\ and N$. In order to solve the dilemma, we have followed 
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L 


iVs 




Ni 


5 


40 000 


200 000 


200 


10 


40 000 


500 000 


500 


20 


6 920 


6 553 600 


400 



TABLE I. Statistics used. For each sample, we run JVhb Heat Bath sweeps and perform Ni measures. Each 
sample is previously thermalized with TVhb iterations. In the L — 20 lattice, we carry out a parallel tempering step 
after each measure. 



the same procedure as in Ref. p4fl to eliminate the bias. One first evaluate the non-linear function with 
the full set of data, then divides the data in two sets of length N\/2 for each of which the function is 
evaluated and the two results are averaged, and the procedure is repeated dividing the data in four sets 
of N\/A measures. We thus have three estimates of the non-linear function, with bias of order 1/Nj, 2/N\ 
and 4/iVi respectively. The three estimates are used in a quadratic (in 1/Nj) extrapolation to l/N\ = 0, 
which is later on averaged over the disorder. In order to have meaningful results from this extrapolation, 
it is crucial that 2r/N\ be a reasonably small number. The values of Nj and Ng for our simulations are 
shown in table |. We remark that we need also to balance the Heat-bath steps, done by the dedicated 
machine SUE, and the parallel tempering steps done by the PC which handles SUE (during this time SUE 
is stopped). 

With our simulation strategy (Ni <C Ns), it is crucial to check that the system is sufficiently thermalized 
while taking measures. A very efficient algorithm for thcrmalizing spin-glass systems is parallel temper- 
ing |l6|-|is||. In order to obtain an efficient parallel tempering we must select a range of (3 values, and 
the number of intervals in this range. The range is fixed in the following way: The faster decorrelation 
time is at the lowest /?; as we run a fixed number of iteration between parallel tempering sweeps, this 
number must be greater than the autocorrelation time at this f3 value. For these values of /3, away from 
the transition point, the correlation length £ is (almost) L-independent. Running around 10 sweeps and 
considering that the correlation time grows as £ 7 [^3j we use finally /3 m i n = 0.70. The largest value must be 
a bit over the crossing point, which we had estimated previously around 0.88. We use then /3 max = 0.92. 
The number of (3 values is fixed by controlling that the probability of changes in the parallel tempering is 
significant. This number depends on L and for L = 20 we have used 12 values of fa obtaining a probability 
of transition around 30%. All the systems are a significant time in the lower (3 values, where decorrelation 
is faster. 

A first thermalization check is summarized in Fig. [j]. The measures taken on a sample are divided in 
twenty blocks, and the correlation-length and the spin-glass susceptibility are calculated with these blocks. 
No thermalization bias can be resolved after the fifth block at the lowest temperature. However, the first 
ten blocks have been discarded for safety. 

Yet the results in Fig. [l] do not really show that we are collecting a reasonable number of measures 
on each sample (so that 2t/N\ is small), because the time needed to obtain a thermalized measure is 
not straightforwardly linked to the time needed to obtain an independent measure. In fact, the latter is 
related with the time needed to overcome free-energy barriers, in order to visit different relevant regions 
of phase-space Q , while the former is related with the time needed to reach at least one relevant region 
of phase space, because the huge number of samples avoids a biased estimation. A better test comes from 
the smoothness of the bias-corrected reweighting extrapolation. In Fig. |^ we show the extrapolated values 
£/L, gi and Xq horn each of the fa of the parallel tempering (alternatively in dotted and dashed lines). 
The mismatch between different extrapolations is much smaller than the error bars, which is due to the 
fact that the same samples are being simulated for all fa values. In fact, for Ni — > oo the mismatch would 
completely disappear, while the statistical error bars would not be smaller unless Ns grows (see Eq. (p2[)). 

Finally, along the simulation the Schwinger-Dyson equations have been checked finding a perfect agree- 
ment. For instance, the equation in |33| holds within a 0.04% for the L = 20 lattice. 

IV. NUMERICAL RESULTS 

Once a set of measures of the finite-lattice correlation-length is at our disposal, the first question one 
can answer regards the nature of the spin-glass phase transition in three dimensions. Indeed, if the 
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FIG. 1. Mean values of the susceptibility and correlation length in the L = 20 lattice for successive bins of 40 
measures, for several f3i values in L = 20 (from top to bottom they correspond to /3u, /3g, Pi and 0o, respectively). 
We remark that /3n = /3 ma x = 0.92 corresponds with our coldest temperature. In the following, we discard the first 
10 bins and average the remaining 10. 
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FIG. 2. Correlation length, susceptibility and <?4 cumulant for the L = 20 lattice, as obtained from the extrapo- 
lation method of Eq. (|l4]) , corrected for bias. We show the results for all the twelve (3 values of the simulation (point 
with error bars) together with the extrapolation (alternatively shown in dashed and dotted lines). The curves have 
been linearly scaled to fit in the figure. 

Kosterlitz-Thouless scenario was realized, (/X on the L — > oo limit, would be zero for temperatures higher 
than the critical one, and then it would abruptly jump at the critical temperature to a finite value. For 
lower temperatures, since the system would still be critical, it would keep a finite value, presumably 
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FIG. 3. Correlation length in units of the lattice size for the 3D Edwards- Anderson model (left), and for the 2D 
XY-model that displays a Kosterlitz-Thouless phase transition. 



growing with lowering temperatures. On a finite lattice, £/L would be a decreasing function of L on the 
paramagnetic phase, and in the critical region (from the critical temperature to lower ones), it would have 
a L independent value, according to the FSS Ansatz, up to scaling corrections. Therefore, in the most 
economic scenario, where scaling corrections are small, the different £/L curves do not cross, but simply 
merge in the low-temperature region. 

In Fig. |§|, we plot £/L for the Edwards- Anderson model in three dimension and for the XY model in 
two dimensions (whose simulation is almost costless in computer time p^j). We see that while the XY 
model follows quite closely the above sketched behavior, the Edwards- Anderson model has a very neat 
crossing. Therefore, one may conclude that the Kosterlitz-Thouless scenario is ruled out by the data, 
unless scaling-corrections of a very exotic nature would be present. Notice, that for the two dimensional 
XY model, £/L and the Binder cumulant behave in the same way with L and T jlO| ]. 

This is an interesting point to compare the behavior of the (74 and G cumulants, defined above, with £/L. 
From Fig. ||, it is clear that the measures for G are much noisier than for £/L. Moreover, also the scaling 
corrections are larger, as made evident by the large shift between the crossing of the 5 and 20 lattice, and 
the crossing of the 10 and 20 lattice. The scaling corrections for G and 34 are of opposite signs, so that 
one can safely conclude that the real critical point is bracketed by this two sets of crossing-points. 

For the critical exponents, our results found using Eq. are displayed in table |l[ Finite-size scaling 
corrections cannot be resolved within errors, specially if one realizes that the results for the (5, 10) pair are 
anticorrelated with the results for the (10, 20) pair (the measure in L — 10 appear once in the numerator 
and the other time in the denominator of Eq.(18)). We take as our final estimate our results for the 
(10, 20) pair, that can be compared with the most recent experimental measures and with other numerical 



calculations displayed in table [II. It would be very interesting to check that systematic errors due to 
finite-size effects are smaller than the statistical ones. This would require to parametrize the scaling 
corrections, and therefore to have precise measures on a wide range of lattice sizes. Unfortunately, the huge 
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T c 


V 


V 


(5,10) 


1.134(9) 


2.39(7) 


-0.353(9) 


(10,20) 


1.138(10) 


2.15(15) 


-0.337(15) 



TABLE II. Critical exponents computed from the crossing points of £/L for (L, IV) pairs. 
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dynamical critical exponent of the available simulation algorithms for this model makes extremely difficult 
to thermalize, at the critical point, lattices larger than L = 20 with present technological capabilities. 

Nevertheless, it is possible to obtain information about u> studying the quotients of cumulants Q g at the 
point where £/ ' L crosses. At this point, Q g should be 1 up to scaling corrections and can be parameterized 
as 



Q g = l + A g L-" . (23) 

Before presenting our results for oj, let us recall the values obtained by Palassini and Caracciolo j[2| 
working in the thermodynamical limit. They computed two different scaling-corrections exponents, A 
and 6. These exponents are obtained form the asymptotic formulae that hold in the scaling region in the 
thermodynamic limit, in the paramagnetic side 

x = A £ 2 -"(l + G(r A )) , £ = B \t\-»{l + 0(\tf)) . (24) 

One can check from the above expressions, that A = u> and 9 = ljv. From their values of A, 8 and v one 
readily obtains 

W (A)=1.3l^ , 0,(0,*/) = 0.78^ . (25) 



We have fitted our measures for g± and g 2 following (|23j). We have 4 points for adjusting 3 parameters, 
but we obtain the value 

uj = 0.84 to.tr > ( 26 ) 

with x 2 — 0.6. See Fig. |5| for more details. Although the statistical error is rather large, the compatibility 
of our result with those displayed in Eq.(25) is reassuring given the difference on the methods. It should be 
also mentioned that the so called analytic scaling corrections have been neglected, which is, a posteriori, 
seen to be a reasonable procedure |39) . 

Finally, we compare our estimate for the critical temperature (T c = 1.138(10)) with that of Ref. |L2| 
(T c = 1.156(15)), the agreement being very good. 
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FIG. 5. Quotients of the cumulants <?2 and (?4 for pairs (L,2L) as functions of (lower figure). The correc- 
tions-to-scaling exponent has been obtained minimizing the \ 2 function (see (|5^)) using the full covariance matrix. 
The horizontal dotted line (in the upper figure) is given by the value of \ 2 at the minimum plus one. 



V. DISCUSSION AND CONCLUSIONS 



We have obtained precise measures of the r\ and v exponents. Moreover, we have done a study of the 
corrections to scaling, in particular we have computed the value of the corrections-to-scaling exponent lu 
in a good agreement with the value reported in Ref. |12|. We remark that our statistical error for the 
77 exponent is 5 times smaller than the experimental error for this exponent and 3 times less than the 
smallest statistical error found in the literature. The fact that our estimates for the (5, 10) lattices and 
the (10, 20) pair coincide within statistical errors gives us some confidence on the smallness of the finite 
size effects, although larger lattices will be needed to make sure that the systematic errors are as small as 
the statistical ones. 

Our comparison with the most recent experimental data is good. The difference between the v 
exponent measured in experiment and our reported value is 0.046(21), roughly two standard deviations. 
We note that the v exponent from numerical simulations is systematically above experimental data. The 
difference for the 77 exponent is 0.023(71). 

The clear crossing of the £/L curves, for different lattice sizes, supports heavily a finite temperature 
second order phase transition and excludes a Kosterlitz-Thouless like scenario (a phase transition of infinite 
order or equivalcntly a line of critical points below the critical temperature). 

Our results for the critical exponents for binary J's agree well with that obtained simulating Gaussian 
couplings H . However the statistical error are still very large and we lack control on the scaling corrections 
(completely in the Gaussian "side"). It might be useful to study of the Gaussian EA model with the 
methodology of this paper but unfortunately the SUE machine is not designed to simulate Gaussian 
coupling. 
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